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Abstract. The isoscalar twist-two generalized parton distributions (GPDs) of the pion and the kaon are 
calculated in a Poincare covariant Bethe-Salpeter constituent quark model. Results are presented for several 
values of the parameters £ and i. The results satisfy the form factor constraints and the polynomiality 
condition. For the pion GPD, also the isospin symmetry constraint is fulfilled. The influence of kinematical 
variables and model parameters on the support of the GPDs is investigated. To this end, the strength 
parameters and quark masses of the constituent quark model are artificially varied. 

PACS. If. 10. St Bound and unstable states; Bethe-Salpeter equations - 12.39. Ki Relativistic quark model 
- 13.60.Fz Elastic and Compton scattering - 14.40.Aq n, K and r\ mesons 



1 Introduction 

Unraveling the substructure of hadrons is a challenging 
issue. Whereas the theory of QCD provides us with the 
equations governing the quark-gluon dynamics, the ex- 
act way in which they form colorless bound states in the 
non-perturbative regime remains elusive to date. It is ex- 
pected that the endeavor of measuring and computing 
generalized parton distribution functions (GPDs) will con- 
tribute significantly to a full description of hadron struc- 
ture pQ. GPDs are a natural unification of form factors 
(FFs) and parton distribution functions (PDFs) within 
one framework. They describe the non-perturbative part 
of the deeply virtual Compton scattering (DVCS) and 
hard exclusive meson production (HEMP) processes. 

A calculation of GPDs from first QCD principles is 
beyond reach at this moment. Therefore, numerous model 
calculations have been presented in the last years [2,3, 
l4ll5l[6ll7l[8ll9l[T0] . In this work, we present results for the 
pion and kaon GPDs computed within the framework of 
the Bethe-Salpeter model developed in refs. [PTII12U13U14] . 
The pseudoscalar ground-state mesons have a simple valence 
quark substructure and possess only one helicity-conserving 
GPD, which makes the calculation less cumbersome than 
for nucleons. Yet, such a calculation provides insight in 
the dependence of the GPDs on hadron structure in the 
non-perturbative regime. 

In order to compare the GPDs computed in a phe- 
nomcnological model with the data from deeply inelas- 
tic scattering experiments, a (^-evolution needs to be 
performed [15]. The evolution equations depend highly 
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on the kinematic region of the process. More specifically, 
one can distinguish between the DGLAP and ERBL re- 
gions, where the Q 2 -evolution is governed by the DGLAP 
and the ERBL equations, respectively. Parton model con- 
straints make GPDs vanish outside these two regions [TBI 
IP?] . A model which resolves the ERBL and the DGLAP 
regions and has a vanishing GPD otherwise, is said to have 
the correct support. 

In ref. [18] we have demonstrated that a support prob- 
lem may arise in dynamic quark models based on the BS 
approach if the interaction depends on the minus compo- 
nent of the relative momentum between the constituent 
quark and antiquark. Note that this condition is necessary 
but not sufficient and that a support problem does not 
necessarily arise when the interaction depends on the mi- 
nus component of the relative momentum. Furthermore, 
we provided representative results computed within the 
explicitly Poincare covariant constituent quark model de- 
veloped by the Bonn group [11] , where a support problem 
did arise. 

In this work, we focus on the influence of the meson 
binding energy on the support. The pion and the kaon 
are ideal for such a study because the former is a very 
deeply and the latter a moderately deeply bound state 
of quark and antiquark. Both the absolute binding en- 
ergy and meson mass are altered by changing the input 
quark masses and effective interaction strength. Yet the 
degree of support violation depends mostly on the rela- 
tive binding energy. Results will be shown for the isoscalar 
twist-two quark GPD of the pion and the kaon at differ- 
ent relative binding energies and for different values of 
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the squared fourmomentum transfer t and relative plus- 
momentum transfer £. 

This paper is organized as follows. Section [2] presents 
the definition of the GPD and the relevant kinematics. 
The Bonn model is introduced in sect. [3J while sect. d]is 
devoted to the calculation of GPDs in this model. In this 
section, we will also introduce the three model variants 
which will be used in sect. [5] There, we show our results 
for the pion GPD. A summary and an outlook to future 
research are given in sect. [6l 



2 Definition and Conventions 




Fig. 1. The Feynman diagram of DVCS on the light cone, with 
the definition of the kinematical variables x and £. 



GPDs are non-diagonal matrix elements of a bilocal field 
operator on the light cone. For partons with spin 1/2 in 
a pseudoscalar meson, the GPD associated with helicity 
conserving partons is defined as follows [16) : 



i r dz~ , 

e 

2tt 



( P'|^(--) 7 +^(f)|P) 



!+=0,: 



(1) 



In this equation, / refers to the flavor of the probed par- 
ton, while P = S^S- and t = A 2 = (P' - P f . Defi- 
nition |lj uses light-cone coordinates. A fourmomentum 
P = if ' i P 1 ) P 2 ) P 3 ) can be written in light-cone coordi- 
nates as p = [p + ,p~ ,p±] with p^ — (p° ±p 3 )/y/2 and 
p ± = (p 1 ,^ 2 ). The skewedness £ and the average plus- 
momentum fraction of the struck parton x are defined in 
fig. [U x denotes the fraction of the average meson plus- 
momentum that is reabsorbed by the meson, while £ is a 
measure for the plus-momentum that is lost in the process: 



P+ - P'+ 
P+ + P'+ 



(2) 



These definitions coincide with the symmetrical variables 
that were introduced by Ji [19] . Their interpretation as 
plus-momentum fractions holds in the infinite momen- 
tum frame. In literature, another set of variables (X, £) 
has been introduced by Radyushkin [5D], with X = |^|| 

and ( = ■ Although Radyushkin's variables are more 
closely related to the ones used in forward kinematics, we 
will use Ji's choice of variables, as they reflect the symme- 
try between the incoming and outgoing hadron states. 

Definition (Q]) is valid when the partons do not transfer 
helicity. It only holds in a coordinate system where q and 
P are collinear and in the three-direction [21] . From defi- 
nition fl} , it is clear that the GPDs are Lorentz invariant 
quantities. This means that, although they are defined on 
the light cone, the GPDs can be calculated in any conve- 
nient reference frame, provided that the average transverse 
hadron momentum Pj_ = 0. Whenever a specific choice is 
needed, we will choose a Breit frame which matches these 
conditions. 



The skewedness £ is restricted to the interval [0, £, m ax) 
with 



£,max 



-t 



AM 2 - t ' 



(3) 



where M is the hadron mass. 

Stringent tests for any model calculation are found in 
the following GPD properties. The first is the relation be- 
tween the electromagnetic form factor and the GPD, 



J dxHf(x,£,t) = F f (t), 



(4) 



where £ factors out and the partial electromagnetic form 
factors are defined through 



(5) 



with e/ the electric charge of a parton with flavor /. The 
second model constraint is called the polynomiality con- 
dition, which states that the nth Mellin moment of the 
GPD is a polynomial in £ of order < n [2T] : 



~ n 

J dx x n - 1 Hf(x ) £,t) = Y,<H(]k)?, 



(6) 



where the coefficients di of the polynomial depend on t. 
The polynomiality condition can be regarded as a more 
general form of the above form factor relation. Finally, for 
pions, isospin invariance and charge conjugation lead to 
the following relation: 

H*(x ) t,t) = -He(-x,Z,t). (7) 

We will come back to these model constraints in sect. [5] 



3 Formalism 

In the Bonn model, mesons are described as bound states 
of a constituent quark and antiquark. The model, based 
on the Bethe-Salpeter equation, is Poincare covariant by 
construction. Explicit covariance is important in the cal- 
culation of the GPD of pseudoscalar mesons for two rea- 
sons. First, the pion is a "deeply-bound" boson whose 
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X 



Fig. 2. Diagrammatic representation of the Bethe-Salpeter 
equation 



static and dynamic properties are best reproduced with 
a relativistic model [22] . Second, in the calculation of dy- 
namic quantities such as the GPD, recoil fourmomenta 
have to be treated relativistically. The Bonn model is a 
few-parameter model with only seven parameters which 
are fitted to the meson mass spectrum [11]. The fact that 
no relativistic corrections need to be implemented is an 
enormous asset to minimize the amount of parameters and 
maximize the predictive power. In refs. [TT ] fT2 ] fT3 ] ll4]. the 
Bonn model is described in detail. In sects. 13.11 and 13.21 
we give a brief summary to make this paper more selfcon- 
tained. 



the constituent quark and the constituent antiquark in 
the meson. Within the instantaneous approximation, all 
retardation effects are neglected. Strictly speaking, this 
is justified for a part of the interaction (e.g. the (colour- 
) Coulomb-interaction) in a specific gauge (e.g. the Coulomb 
gauge) only. Otherwise this approximation will lead to a 
non-local and non-causal model. Nevertheless, the model 
equations lead to frame-independent results. Hence, the 
model can be called Poincare covariant. The second as- 
sumption is that the full quark propagators Sf (p) can 
suitably be approximated by free fermion propagators with 
an effective constituent quark mass rrij, 



Sf( P ) ^ 



(12) 



The constituent quark masses are model parameters. Be- 
cause isospin symmetry is exact in the Bonn model, the 
calculation of the low-lying meson spectrum requires only 
two mass parameters: the non-strange and the strange 
quark mass. 



3.1 The model ingredients 

A meson with on-shell fourmomentum P is described by 
the Bethe-Salpeter amplitude 



(8) 



where T is the time ordering operator acting on the Heisen- 
berg fermion field operators ip a - The indices in Dirac, 
flavor and color space are combined in the multi-indices 
a = (a, /, c). The variables x\ and xi are the fourvectors 
denoting the spacetime positions of the quark and the an- 
tiquark, respectively. 

The Bethe-Salpeter amplitude (O is the solution of the 
Bethe-Salpeter equation, which in momentum space reads 
[H[23]: 

XP = -iG pKpxp ■ (9) 

Here, indices and arguments have been suppressed for no- 
tational simplicity. It is tacitly assumed that one inte- 
grates over arguments and sums over indices that occur 
twice. G p is the product of the one-particle propagators, 
while Kp denotes the interaction kernel. The pictorial rep- 
resentation of eq. ([9]) is shown in fig. [2] The normalization 
of the Bethe-Salpeter amplitudes is given by [T2"] : 



XP 



d 



dP^ 



(G-^ + iK P ) 



Xp = 2iM 2 



(10) 



p=p 



Note that eq. (|10p is written in a frame independent way. 

To transform the Bethe-Salpeter equation © into a 
solvable integral equation, two Ansatze are made. The first 
assumption is the instantaneous approximation of the in- 
teraction kernel Kt 



K P (p,p') = V(p ± p,p' J _p), 



(11) 



where p±p = p — (p ■ P/P 2 )P is the fourvector perpen- 
dicular to P with p the relative fourmomentum between 



3.2 Reduction to the Salpeter equation 

The Salpeter amplitude <P(p) in the rest frame of the 
bound state is defined as follows: 



f dp 

®(p) = / -^rXp{p°,p)\p=(M,o) 



(13) 



Making use of eqs. (fTT]) and (| 13[) . the p°-integration in eq. 
© can be carried out, leading to the well-known Salpeter 
equation [2"1] : 



Hp) 



d V A[ (p) 7 ° [V[p, p')$(p')]^Aj (-p) 

(2tt) 3 M + UJl + U! 2 

dV A+(p) 1 ° [V(p, p')$(p')h°A 2 (-p) 
(2tt) 3 M - u>i - w 2 



(14) 



with the energy projection operators = (uij±Hj)/(2ujj), 



the Dirac Hamiltonian Hj(p) 
energy uj 



7° (7 ' P + 171 j) an d the 

ij + \p\ 2 for the jth constituent quark. 

The potentials used in our calculations are the con- 
finement interaction V and the 't Hooft instanton induced 
interaction Vm [25] , The confinement potential rises lin- 
early with the interquark distance and is multiplied by a 
Dirac structure which reproduces the observed mass split- 
tings in the meson spectrum [26U27j : 



V(|x g - Xq\) = i(a c + b c \x q - Xq\) 

x (1 ® 1 - 7 5 ® 7 5 - 7m ® 7^) , (15) 

where x q ^ is the position of the constituent (anti-)quark. 
The instanton interaction accounts for the mass splittings 
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in the pseudoscalar and the scalar sectors [22] 
turn space, it can be written as 

/ J^V I u{j P ,p'W)=±G{g,g') 
dV 



In momen- 



(2tt) 3 



n A {p 1 p')(lTr[<P(p l )]+ 1 5 Tr[ 1 5 <P(p')]) 



(16) 



where G(g,g') includes the flavor dependent couplings g 
and g' , and TZa is a Gaussian regulating function with cut- 
off A [26II27] . The 't Hooft instanton induced interaction 
can account for the low mass of the pion. The confine- 
ment and the instanton interaction contain five additional 
model parameters, bringing the total to seven. These seven 
model parameters are fitted to the Regge trajectories and 
the pseudoscalar ground state masses and are kept fixed 
in the calculation of dynamic observables. In this sense, 
the results in sect. [5] are predictions. 

The Salpeter amplitudes <P are calculated by solving 
eq. (|14|) . Imposing the normalization condition, written 
in terms of the Salpeter amplitudes, then yields the mass 
spectrum jl"2] . 



4 GPDs in the Bethe-Salpeter Quark Model 

The first step in the calculation of the generalized par- 
ton distribution defined by eq. (Q]) is the calculation of 
the bilocal current matrix element {P'\ip(— §)7 +- 0(f )|-P)- 
In a Bethe-Salpeter based approach, the matrix element 
of any dynamic variable can be calculated with the Man- 
delstam formalism [28]. This formalism acts as a starting 
point for the derivation of the GPD in terms of the Bethe- 
Salpeter amplitudes. A straightforward calculation yields 
the bilocal current matrix element in lowest order 1271: 



= a 



A «Tr jxp,(4, |) ((« @x' 2 -mi)xp(4)-|)) 7 + } 

+ J d 4 yi zTr{((z fa - TO2 )xp,(-| )2 /0) 7 + Xp(|,yi)} • 

(17) 

The second term in eq. (|17[) (coupling to the quark) is 
depicted in fig. [3] Inserting eq. (fl~T|) in eq. (Q]) gives 

H(x, C, t) = m(x, C, t) + H*(x, £, t) , (18) 

with the term originating from the quark-current coupling, 



El 



v 




Fig. 3. Schematic representation of the second term in expres- 
sion (|17p for the plus-component of the bilocal current matrix 
element. 



and the term originating from the antiquark-current cou- 
pling, 



P' 



(20) 



Equations (fT71) - ([2U|) make use of the full Bethe-Salpeter 
amplitude \- This amplitude can be reconstructed from 
the Salpeter amplitudes as follows: once the Salpeter equa- 
tion is solved, the vertex functions Tp — G~pXp can be 
calculated. In the meson rest frame, these vertex functions 
read: 

rp(p xP )\p HMt0) = r( P ) = -iJ ^JL[v{p, P 'W)] . 

(21) 

By Lorentz boosting the resulting Bethe-Salpeter ampli- 
tude from the meson rest frame to the frame in which the 
meson has on-shell momentum P, the amplitude in the 
latter frame is found: 



Xp(p) = Sa p X(m,o)( A p 1 P) S aI 



(22) 



In this equation, Ap denotes the Lorentz transformation 
and Sa p denotes the corresponding boost operator acting 
on the fermion field operators in eq. ([9]). Due to the instan- 
taneous approximation, the boost properties of the Bethe- 
Salpeter amplitude are described by the above equation. 
All interaction dependence effectively enters the boost op- 
erator through the bound state's mass. 

Written in terms of these vertex functions, the quark 
GPD arising from the Mandelstam formalism contains 
three quark propagators Sf- 



2 J (2tt)4% 2(1+0 P 



Tr {rp,(p+^)S F (^+p+jh + 



(19) 



(23) 

These propagators can be directly linked with the inter- 
mediate quark lines in the diagram of fig. [3] It turns out 
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Parameter 


Full model 


Reduced model 


IQM model 


m n [MeV J 
m s [MeV] 


o on 
380 

550 


o on 
380 

550 


800 


dc [MeV J 
b c [MeV/fm] 


1 1 or 

-lloo 
1300 


line 

-lloo 
1300 


1 1 9 K 

1300 


y l>*e v j 
g' [GeV- 2 ] 
A [fm] 


1.62 
1.35 
0.42 


0.0 
0.0 


1.62 
1.35 
0.42 


Mtt [MeV] 
[MeV] 


141 
619 
4.39 


572 
188 
0.33 


1095 
505 
0.46 


Mk [MeV] 
[MeV] 
A M /M K 


506 
424 
0.84 


728 
202 
0.28 





Table 1. Overview of the parameters of the three models that 
were used in this work: the constituent quark masses, the con- 
finement offset and slope, the 't Hooft interaction range and 
the 't Hooft interaction strengths. Also presented is a summary 
of the masses M, binding energies Am and relative binding en- 
ergies Am /M of the pion and the kaon in the different models. 



that the denominators of these three propagator terms en- 
sure the correct support region for the GPDs, x G [— £, 1] 
(ERBL and DGLAP regions) for p~ independent vertex 
functions. For p~ dependent vertex functions, it is a priori 
unclear whether the GPD will be confined to the support 
region. In previous work, we have shown that the Bonn 
model is prone to a support problem 18] . In the next para- 
graphs, we will quantitatively investigate which physical 
parameters influence the support behavior of the model. 
The role of the variables £ and t will be analyzed, and 
the dependence on the binding strength will be examined 
through a comparison of the GPD of the kaon and the 
pion in three different model variants. 



"deeply bound" . The kaon is moderately bound in this 
model, with a mass of M£ u " = 506 MeV and a mass de- 
fect of A f M* = (m„ + m s - M f K uU ) = 424 MeV. 

4.1.2 Reduced model 

In the reduced model, the 't Hooft instanton induced in- 
teraction is omitted. As we have mentioned in sect. 3.2, 
the instanton interaction accounts for the deep binding 
of the pion (and to a lesser degree also of the kaon) . Ac- 
cordingly, neglecting the 't Hooft interaction will provide 
insight into the importance of binding effects in the GPD 
results. As a matter of fact, the calculated pion mass in- 
creases to Ml ed = 572 MeV in this approach (A M d = 
188 MeV). The kaon mass increases to M^ d = 728 MeV 
{A'* = 202 MeV). 



4.1.3 Increased quark mass (IQM) model 

Not only the 't Hooft instanton induced interaction has 
an effect on the pion binding energy. Also the non-strange 
constituent quark mass m n affects the pion mass and mass 
defect. To investigate the influence of the binding energy 
on the support of the generalized parton distributions, 
both mechanisms must be studied. 

In the IQM model, we will only show results for the 
pion GPD. After combining the IQM model results with 
the kaon results from the full model, one can determine the 
influence of the heavy quarks. Increasing the non-strange 
quark mass by more than a factor of 2 to m n = 800 MeV 
yields a pion mass of M^ M = 1095 MeV. The mass defect 
in this model is A[^ M = 505 MeV. 



4.1.4 Model summary 



4.1 Model variants The masses, binding energies and relative binding energies 

(defined as the binding energy devided by the mass) cal- 
In the forthcoming section, the GPDs of pseudoscalar mesons culated in the different models are summarized in table 



will be shown in three different model variants: the full 
model, the reduced model and the increased quark mass 
(IQM) model. The parameters of these models are pre- 
sented in table [TJ Notice that the parameters a c and b c 
of the confinement interaction remain fixed for all three 
model variants. 



4.1.1 Full model 

The full model is the one referred to as Model B in ref. 
|26j . It provides an accurate description of the pion and 
other meson properties such as its mass, electromagnetic 
form factor, electroweak decay widths, etc. In the full 
model, the pion mass is calculated as Af/ U " = 141 MeV 
which implies a large mass defect of Am" 1 = (2m n — 



M ful, 



619 MeV; accordingly we shall call the pion 



[U The deep binding of the pion in the full model is re- 
flected in the high relative binding energy. Am /M in the 
full model is much smaller for the kaon than for the pion. 
Further, Am/M of the kaon in the full model is larger 
than in the reduced model, and also larger than Am/M 
of the pion in the reduced and the IQM models. In sect. 
5, we will come back to these (relative) binding energies. 



4.2 Model constraints 

In sect. 2, we introduced three constraints that serve as 
stringent tests for any GPD calculation. These are the 
isospin symmetry relation ([7]) for the pion GPD, the form 
factor relation (j4j and the polynomiality condition (JH) . In 
this section, we elaborate on these constraints, and show 
that they are fulfilled in our model. 
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The pion up and down quark GPDs must fulfill the 
isospin symmetry relation of eq. I|7|. The equality is ex- 
act in our calculations. The strange quark content of the 
kaon prevents an isospin symmetry relation of the type 
([7]). written. The results show a small difference between 
the quark and antiquark GPDs at opposite x {e.g. in the 
full model, H£ + (x = 0.5, £ = 0, t = -0.5 GeV 2 ) = 0.764 
versus H s K+ {x = -0.5, £ = 0,t = -0.5 GeV 2 ) = -0.784). 
We will elaborate on these differences in sect. [5] 

The GPD is related to the electromagnetic form factor 
through relation Taking into account eq. (O, which 
relates the partial form factors with the meson form factor, 
and the isospin symmetry relation |J7_J) for the pion, one 
finds that 



dxH2 + (x,S ) t)=E r+ (t). 



(24) 



Note that, due to the support properties of the GPDs in 
the Bonn model, the integration domain is x € (-co, +oo). 
For the kaon, the relation becomes 



/+oo 
dxH%-+(x,£,t) 
-OO 

/+oo 
dxH s K+ (x,£,t)=F K+ (t). 
-OO 



(25) 



where e u — 2/3 and e s = —1/3. 

We have compared the results of eqs. (|2"4"|) and 
with a direct computation of the electromagnetic form 
factors i 7 ' 7r + and F K + in the Bonn model [11] . These nu- 
merical calculations were performed independently and 
yielded results which were compatible at the few-% level 
[27] . When discussing these form factors results, it is worth 
stressing that the instantaneous approximation does not 
result in current non-conservation. Although the EM cur- 
rent in lowest order in the instantaneous approximation 
generally does not fulfill the Ward-Takahashi identity, it 
does so in the case of a 0" — * 0~ elastic matrix ele- 
ment [14] . 

An even more stringent test than the form-factor con- 
dition of eq. (j4|) , is the polynomiality condition of eq. ([6]) . 
This condition was verified numerically for different values 
of t up to order n = 5 [27] . 

In [29], the authors point out that the parton distri- 
bution (PD) of the pion should fulfill four conditions: (a) 
the support region is — 1 < x < 1, (b) the first moment of 
the GPD must be 1 (normalization) for any value of /xi, 
(c) in the case of the pion the second moment of the PD 
must be 1/2 (momentum conservation) and (d) isospin 
and momentum conservation imply that the PD of the 
pion must be symmetric under the change x — ► (1 — x). 
As to condition (a), the pion vertex function in our model 
is expressed in Minkowski spacetime and has an intricate 
analytic structure as opposed to Ref. [29] . It was already 
argued in our previous work [18] that this can lead to a 
support problem. Condition (b) is theoretically fulfilled in 
our model and has been checked numerically [27]. Condi- 
tion (c) is not fulfilled in our model and may be attributed 



to large off-shell effects. Concerning condition (d), isospin 
symmetry translates into condition ([7]) in our model, and 
this has been confirmed numerically to be valid for every 
value of £ and t. 



5 Results and Discussion 

The 7r + up-quark GPD results in the three different model 
variants are shown in figs. H] - [U the K + up-quark and 
strange-antiquark GPDs in figs.[7]-[8l A three-dimensional 
picture of the (full model) ir + up-quark GPD as a function 
of x and t with skewedness £ — is presented in fig. [9] 
For all cases, the results are shown for a representative 
selection of t and £ values. 

It was shown in ref. [18] that GPDs in the Bonn Model 
violate the support condition. This can also be deducted 
from figs. S]-[5] In general, the curves tend to have longer 
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Fig. 4. The pion GPD H™+ in the full model for different 
values of t. Values for £ shown are £ = (red, solid line), 
£ = 0.4 (green, dashed line) and £ = 0.7 (blue, dotted line). 



tails as \t\ decreases. Especially for the deeply bound pion 
in the full model, this effect is clearly visible (see fig.[9j). As 
the relative binding energy decreases, the effect becomes 
smaller. Furthermore, whereas the maximum of the GPD 
curve lies outside the support region for the deeply bound 
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Fig. 5. The pion GPD H™+ in the reduced model for different 
values of t . Values for £ shown are £ = (red, solid line) , £ = 0.4 
(green, dashed line) and £ = 0.7 (blue, dotted line). 



pion in the full model (peak value at \x\ > 1), it lies within 
the support interval for the other pion and kaon calcula- 
tions. 

Another observation which can be made, is that the 
pion GPD in the full model displays a clear shoulder in the 
region x € [0, 1] at \t\ < 0.1 GeV 2 (see fig. gj). This shoul- 
der is much less pronounced or even absent in the other 
results. In this respect, it is interesting to note that the 
computed pion form factor shows a bump at low Q 2 = —t, 
which can be ascribed to the instantaneous approxima- 
tion [11,30,31]. The observed shoulder in the pion GPD 
might therefore be an artefact of this Ansatz. 

A direct comparison with other model calculations ([21 
[4]) is not straigthforward due to the different behavior 
with respect to the support properties. A double-peaking 
behavior for x = (1 ± £)/2 in the reduced or IQM model 
variants is not seen. In both of these models, however, 
the binding energy represents a significant portion of the 
meson mass: about 1/3 (1/2) for the reduced (IQM) model 
(see Table [1]). This is in contrast with the calculations 
of Ref. [5], where the meson mass is very close to the 
sum of the quark masses and the double-peaking behavior 
emerges. 



Fig. 6. The pion GPD H™+ in the IQM model for different 
values of t. Values for £ shown are £ = (red, solid line), 
£ = 0.4 (green, dashed line) and £ = 0.7 (blue, dotted line). 



The observed model dependence of the GPD properties 
points at a significant relative binding energy dependence 
of the support of the generalized parton distributions in 
the Bonn model. To quantify the support problem, a sup- 
port parameter (0) is introduced. For the pion, the gen- 
eralized quark and antiquark distributions are related via 
the isospin symmetry relation |[7J), so that the knowledge 
of one of them implies the knowledge of the other. The 
support parameter is therefore defined through the quark 
GPD: 



ftt\H2 + (x,Z,t)\dx 



(26) 



<fi is a measure for the relative support violation: when 
the correct kinematical regions are resolved, = 1. The 
smaller <fi, the worse the support. For the kaon, the isospin 
symmetry relation |(7j) is not valid and different support 
parameters are introduced for the quark and antiquark 
GPDs: 



A. 

r i 

J — oo 1 



\H» + (x,£,t)\dx 
H* + (x,£,t)\dx 



(27) 
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Fig. 7. The kaon GPDs H™ + and H a K+ in the full model for 
different values of t. Values for £ shown are £ = (red line), 
£ = 0.4 (green line) and £ = 0.7 (blue line). The solid line 
refers to the u GPD, the dashed line to the s GPD. 



Fig. 8. The kaon GPDs and H a K+ in the reduced model 

for different values of t. Values for £ shown are £ = (red line), 
£ = 0.4 (green line) and £ = 0.7 (blue line). The solid line refers 
to the u GPD, the dashed line to the s GPD. 



and 



r i 

J — oo 1 



(28) 



The difference between <j) q and <ffi will be used to investi- 
gate the flavor dependence of the support problem. 

Table [2] lists the values of the support parameter <j> 
for the pion GPD results of figs. [4]- [6] Table [3] shows the 
values of the parameters <p q and <jfl belonging to the kaon 
curves of figs. [7] and [8] 

For the pion GPD, it is seen from table [2] that the 
support is much better in the reduced and the IQM mod- 
els than in the full model. In going from the full model 
to either of the two other models in which the pion is 
less deeply bound, the support parameter increases signif- 
icantly (e.g. from 0.08 in the full model to 0.60 in the re- 
duced and 0.66 in the IQM model for t = -0.01 GeV 2 and 
£ = 0). We observed in the previous paragraph that, espe- 
cially for the pion in the full model, the GPD curves tend 
to broaden with decreasing \t\. Table [2] shows that this be- 
havior can be attributed to the support properties of the 
curves. For fixed skewedness £, the support improves for 
increasing |i| except for very large values (|£| > 5.0 GeV 2 ). 
For a fixed |t|, the <j> exhibits hardly any dependence on £. 



Also for the kaon GPD, the support improves when the 
instanton induced interaction is switched off. Notice that 
the values of the two support parameters <p q and <\fl are 
similar. This result is compatible with the finding that the 
quark and antiquark GPD have similar shapes and shows 
that the small differences between both GPDs do not alter 
the support significantly. Combining these observations 
with the fact that the support does improve significantly in 
the IQM model with respect to the full model, we conclude 
that concerning the support the binding strength is more 
important than the particular constituent quark mass or 
flavor. 



6 Summary and Conclusions 

In this work, the isoscalar twist-two generalized parton 
distributions of the pion and the kaon were calculated in 
the Poincare covariant Bethe-Salpeter constituent quark 
model developed by the Bonn group. The first moment 
of the GPD in the full cc-region equals the electromag- 
netic form factor. Results were shown in the full x-region 
for several values of £ and t. It turns out that the Bonn 
model violates the support condition. We have illustrated 
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Table 2. Values of the pion support parameter <f> of eq. (|26[l corresponding to figs. [4]- [6] 











t (GeV^) 






Parameter Model 


t 


-0.01 


-0.05 


-0.1 


-0.5 


-1.0 


-5.0 


4> q Full model 





0.46 


0.48 


0.50 


0.58 


0.61 


0.47 




0.4 








0.60 


0.62 


0.48 




0.7 










0.56 


0.56 


Reduced model 





0.68 


0.69 


0.70 


0.73 


0.73 


0.64 




0.4 








0.76 


0.75 


0.62 




0.7 












0.63 


<f> q Full model 





0.43 


0.45 


0.47 


0.56 


0.61 


0.42 




0.4 








0.58 


0.61 


0.43 




0.7 










0.55 


0.57 


Reduced model 





0.65 


0.66 


0.67 


0.72 


0.73 


0.65 




0.4 








0.73 


0.74 


0.64 




0.7 












0.69 



Table 3. Values of the kaon quark and antiquark support parameters of eq. (|28[) corresponding to figs. 0- [8] 



the strong correlation between the support and the (rel- 
ative) binding energy of the meson. Therefore, the deep 
binding of the pion induces strong support violations. We 
have found that £ hardly influences the support, whereas 
a moderate dependence on t is predicted. The constituent 
quark masses have an impact on the support, but only 
through the (relative) binding energy. A mass difference 
between constituent quarks hardly affects the correspond- 
ing GPD curves. 

The reason for the violation of the support condition in 
the Bonn constituent quark model is not fully understood: 
It might be related to the instantaneous approximation 
made. Although this assumption is formally covariant, the 
neglect of retardation effects could lead to wrong analytic 
properties of the hadronic vertex functions, in particular 
for "deeply-bound" states such as the pion. Some hints 
in this direction were already found earlier in the cal- 
culated pion electromagnetic form factor, which showed 
a conspicuous unphysical structure at small momentum 
transfers [301I3T]. The present detailed analysis provides 
additional support for this conjecture. 

The observed support properties of the Bonn constituent 
quark model could also be interpreted as follows. As the 
Bonn model is a soft scale model, with neither on-shell 
particles nor asymptotic freedom properties, the support 
violation can be considered as a manifestation of soft scale 
physics. This prohibits the use of the available QCD-based 



evolution equations to the model. In contrast, a soft-scale 
model with the correct support properties can be evolved. 
It should be clear, however, that also in this case, the use 
of the evolution equations relies on an extrapolation of 
perturbative QCD to soft scales. Additionally, it should 
be noted that the Mellin moments of the GPDs at £ = 
can be evolved, even in the case of a support violation. 
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